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(N| i Abstract. Using ID NLTE radiative hydrodynamics we model the influence 

^ . of the particle beams on the Ha line profile treating the beam propagation and 

& 1 the atmosphere evolution self-consistently. We focus on the influence of the 

non-thermal collisional rates and the return current. Based on our results, we 
propose a diagnostic method for determination of the particle beam presence in 
' the formation regions of the Ha line. 

CO 

^ ' 1. Introduction 

o ' 

Some of the flare models assign a fundamental role to the high energy particle 
beams in the flare energy transport. As the beams interact with the ambi- 
ent plasma, their energy is dissipated and transformed mainly into the ther- 
mal energy of the transition region and chromosphere plasma. Several mod- 
els studied electron and proton beams as heating agents (e.g., iMariska et al 



19891: Emslie et al.l [l998l as well as their influence on spectral line profiles (e.g., 



Alfred et alJlSh . The propagation of electron beams is ine vitably connected 



with the so-called return current (RC, Ivan den Oordl Il990l ) which also con- 
tributes to the beam energy dissipation. Besides the heating, the beams in- 
fluence atomic level populations of the ambient plasma via collisional excitation 
and ionization. The work presented here concentrates on these two effects which 
are commonly neglected in the flare modeling and assesses their importance on 
rN ' the formation of Ha line in solar flares. 



2. Model 

The beam propagation and energy deposition is modelled by a test particle 
approach consistently with the hydrodynamics of the atmosphere and NLTE ra- 
diative transfer in the transition region, chromospher e, and photosphere. De tails 
of the model and the methods used are described in Kasparova et al. I (l2005h. 



We study the response of quiet Sun atmosphere (VAL C from lVernazza et al. 
Il98ll ) to beam pulses of short duration, 1 s with sinus-like time modulation, and 
power-law energy spectrum with 5 = 3. 

3. Hydrodynamics and Beam Propagation 

The model takes into account Coulomb collisions of the beam with ambient 
neutrals and electrons, scattering of beam electrons, and optionally RC for the 
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Figure 1. Total energy deposit (solid) and energy deposit to hydrogen 
(dashed) at beam flux maximum. Top: electron beam with and without 
RC included. Bottom: proton beam with E 1 = 5 MeV and E x = 20 MeV. 



case of an electron beam (return curre nt is a factor of E e jE v lower for protons 
carrying the same power as electrons (jBrown et al.lll990l ) and is neglected). The 
return current is included in a runaway approximation assuming a — 0.1, i.e. 
10% of ambient electrons carry RC . For details and other ap proximations of RC 
in solar atmosphere conditions see IVarady et al. f2005. l2007h . 

Figure [1] shows that RC significantly increases the energy deposit of electron 
beams at heights > 1500 km leading to corresponding increase of the tempera- 
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Figure 2. Time evolution of the Ha line center (left column) and wing 
(right column) intensities. Top: electron beam for F = 10 10 erg cm -2 s _1 and 
(5 = 3 with and without RC. Bottom: proton beam for F = 10 12 erg cm -2 s _1 , 
5 = 3, and E\ = 5, 20 MeV. The solid curves refer to cases with C nt included, 
dashed lines to cases without C nt . 



ture and ionization. The atmospheric response to proton beams was modelled 
for two different values of the low-energy cutoff E\ — 5, 20 MeV. Note that deka- 
MeV protons produce approximately the same amount of hard X-rays as deka- 
keV electrons (lEmslie fc Brownlll985l ). Figure [D illustrates that proton beam s 
deposit their energy into deeper layers than electron beams (lEmslie et al.| [l996). 
The energy deposit of proton beams with lower value of E\ peaks at higher atmo- 
spheric layers and is larger at heights above in comparison with a proton beam 
of E\ — 20 MeV. Consequently, the temperature at these layers is substantially 
increased. 



4. Ha Line Profiles 



The hydrogen level populations are affected also by non-thermal collisional ion- 
ization and excitation. The corresponding r ates C nt are dir ectly proportional to 
the b eam's energy depos it on hydrogen (see lFang et al.|[T993l for electron beams, 
and iHenoux et all [T993 for proton beams). Their influence o n the tempora l 
evolution of the Ha pr ofile was studied for electron beams bv iHeinzell dl99lft : 
Kasparova et al.1 (|2005l ) . Here, we describe in detail their effect for proton beams 
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Figure 3. Ha contribution function CF for an electron beam. RC and C ] 
were included. Solid curve: energy deposit on hydrogen. 



and the influence of RC for electron beams. The temporal evolution of the Ha 
line-center ( AA = A) and win g ( AA = 1 A) intensities are shown in Fig. [21 As 
discussed in Heinze 1 (Il99lh and Kasuarova et al. (2005) for the case of electron 
beams, C nt cause a decrease of the line-center intensity at the very start of the 
beam propagation and enhance the wing intensity later on (mainly for fluxes 
J> 10 11 erg cm -2 s _1 ). The line behaviour can be understood in terms of contri- 
bution functions CF to the outgoing intensity given by I\ = f CF\dz, where z 
is the height. 

Figure [3] shows the evolution of Ha CF for the case of an electron beam 
with RC and C nt included. A decrease of the line center intensity is caused by 
an increase of opacity due to C nt (see CF at 0.2 s). Later on, a new region of 
wing formation occurs at the layers where the energy deposit peaks (see at 
0.4 s). Such a layer is not formed if C nt are not considered. Since RC causes 
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Figure 4. Comparison of Ha contribution function CF for proton beams 
with Ei = 5 MeV (first column) and E x = 20 MeV (second column). C nt 
were included. 



446 



Kasparova et al. 



heating of the top parts of the atmosphere, it is responsible for the increase 
of the line center intensity at 0.4 s (see Fig. [2j) forming in a narrow region at 
- 2000 km (see Fig. Ej). 

The situation for a proton beam is shown in Fig. [H In the case of E\ — 
20 MeV, the beam energy is deposited in regions where the line wings are formed. 
If C nt are included, they cause increase of opacity which leads to a drop of both 
the wing and the line-center intensities (see CF at 0.2 s). Later on, a new wing 
region occurs. Neglecting C nt , the line intensity does almost not change because 
the temperature structure is not significantly affected by the beam propagation. 

However, for a proton beam of E\ — 5 MeV, the temperature increase is 
large enough to form a new region of the line center intensity, as in the case of 
RC for an electron beam. Also, the energy deposit is large enough to create a 
new layer of strong wing intensity (see CF at 0.4 and 1.1 s and Fig. ED- The 
second peak of the line center intensity at about 1.1 s is due to broadening of 
the formation region. Similarly to the electron beams, the line wings are not 
changed if C nt are not considered. 

5. Conclusions 

The Ha line is influenced by temperature structure resulting from the beam 
propagation and return current as well as the non-thermal collisional rates. We 
propose to use the decrease of Ha line-center intensity as diagnostic indicating 
the presence of particle beams. We also predict that proton beams with deka- 
MeV low-energy cutoffs produce only decrease of the Ha in comparison with 
quiet-sun intensities. 
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